Peierls substitution in the energy dispersion of a hexagonal lattice 
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The method of the Peierls substitution in studying the magnetic subband structure of a hexagonal 
lattice is re-examined. Several errors in the formalism of a couple of recent papers are pointed out 
and rectified so as to describe the effect of the magnetic field pertinently. 



PACS numbers: 71.28.+d, 71.20.-b, 73.20.Dx, 71.45.Gm 



Recently, Fekete and GumbsB reported on the study 
of the energy spectrum of two-dimensional (2D) Bloch 
electrons subject to a perpendicular sine-modulated mag- 
netic field 



B z {x,y) = B +B 1 {x)=B + B l sm{2nx/T x ) (1) 

by making the Peierls substitution k — > (p + eA)/h in the 
first-order energy dispersion e(k) obtained by the tight- 
binding method. The content of Ref. 1 consists of two 
parts. The first part deals with the formalism that de- 
scribes the effect of magnetic modulation (i.e., B\ ^ 0) 
on the energy spectrum of a square lattice, where a cor- 
rection of some errors in a papera was given. (Concerning 
this problem, it may be worthy of consulting Ref. 3.) The 
second part deals with the energy spectrum of a hexago- 
nal lattice in the presence of magnetic modulation, where 
the authors attempted Jto generalize the formalism pre- 
sented in another paperQ which deals with the problem of 
the energy spectrum of the hexagonal lattice under a uni- 
form magnetic field. Through the studies, the authors of 
Refs. 1,2, and 4 calculated energy eigenvalues, in the ab- 
sence and/or presence of magnetic modulation, by diago- 
nalizing the effective Hamiltonian obtained by the Peierls 
substitution. However, as for the hexagonal lattice, the 
formalism and the numerical data presented in Refs. 1 
and 4 are flawed in several ways, whose origin lies in a 
simple but fundamental mistake made at the early state 
of the Peierls substitution. In this comment, we shall 
point out the mistake and the relevant errors of Refs. 1 
and 4, and derive an exact formalism that enables to 
pertinently describe the effect of the magnetic field given 
by (|l|) on the energy spectrum of the hexagonal lattice. 
For illustration purpose, we will present numerical data 
of the energy eigenvalues obtained by diagonalizing the 
corrected matrix with and without magnetic modulation. 

To derive the Hamiltonian matrix that describes the 
effect of the magnetic field given by ([j]), we consider the 
energy dispersion of a 2D hexagonal lattice given by 



e(k) — 2 {to cos(k x a) + t + cos(k x b + k y c) 
+t_ cos(k x b — k y c)} , 



(2) 



where a is the lattice constant, b — a/2, and c = \^3a/2. 
Making the Peierls substitution in (||) yields the effective 
Hamiltonian 



H = to<Z iv * a l h + t + Q i [P^ b +Pv c + eA v( x ) c ]/ h 

-\-t-c i ^ Pj:b ~' PyC ~ eAy( ' x:ic ^ h + c.c, (3) 



where 



A y (x) = B x - (SiT x /27r) cos(2nx/T x ) 



(4) 



is the y component of the vector potential under the Lan- 
dau gauge. 

The four exponents of the form e Apx+Bpy+f< - x ^ in (||) 
contain the operators p x and x simultaneously in their 
arguments, unlikely to the case of the square lattice, and 
thus we should be very careful in manipulating them. 
Since p y commutes with p x and x, we can easily decou- 
ple e A P*+ B Py+f( x ) as 



2 A Pl + B Py + f{x) _ e B P y c A P:c +f(x) _ c Ap x +f(x) e Bp y 



(5) 



However, since p x does not commute with x, we cannot 
simply decouple e Ap,c+ f( x > as 



e A P * + f{x) _ e A P x e f{x) or e /(x) e Ap x 



Instead, it should be written as 

c Ap x + f(x) = e Ap ae fo f( x + ihAX ) dX 



(6) 



(7) 



which is a generalization of the Weyl's formula, and can 
be derived by using the Baker-Hausdorff lemmad. In Par- 
ticular, for f(x) — Cx + Dcos(Kx), (]?]) becomes 

e Apx + Cx+Dcos{Kx) _ e Ap Xe C(x+ihA/2)-iD£(x)/hAK ^ 

where 

£(x) =sm[K(x + ihA)] -sin(Kx). (9) 
Combining (||) and (|J) with (||) , we can write 

n = t (e ip * a / R + e - ip * a/h ^ 

4-t + c~ iva / 2 (^Q i Py c / n Q i P^ b / Tl Q i (' Rax ' + v -) 
_|_ g— iPyC/ft-Q— iPxb/frg— i(nax— v + ) 

+i_e" ra / 2 ^ e - i Py c / h e i Px b / h e - i ( 7Tax '+ l/ -) 

_j_ e ipyc/n e -ip as b/Ti e i(-Rax : ' -v+)^j ^Q) 
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FIG. 1. Energy eigenvalues versus a for 7 = with (a) 
(to, *+,*-) = (1,1,1) and (b) (t ,t+,t_) = (2,1,1). Cal- 
culations are performed for a = p/101 (1 < p < 100) and 
(*»,*») = (o,o). 

where 

f±(a0 = {sin[7r/3(a;' ± 1)] - sm(n0x')} , (11) 
7rp z 

and 

a = eB ac/2Trh, (3 = a/T x , 7 = B 1 /B , x' = x/b, (12) 

Denoting the lattice points as (m, n) — (x/b,y/c) and 
using the translational property e~ l ^'^\r) = \f + £), we 
can write the Schrodinger equation (r\H\ip) = Eijj(r) as 

Elpm,n — to {lpm-2,n + ^m+2.n) 

+t + (e^^-^m-l.n-l +e lBm ^m+l^+l) 

+t_ (e ie — 1 ^ m _i, n+ i + e- i( ^ m+1> „_i) , (13) 

where 

m (= 0<? + 0&>) =7ra(m + l/2) 

— ^ {sin[(m + 1)tt/3] - sin(m7r/3)} (14) 
7r/!r 

is the magnetic phase factor. Note that 9 m plays a key 
role in studying the magnetic subband structure within 




-8 1 1 1 1 1 

1 
p/q 



FIG. 2. The same as figure 1, except that 7 = 2 and T x = 4. 

the tight-binding approximation since all the effects of 
the applied magnetic field are held in it; 9m reflects the 
effect of Bq while 8m reflects the effect of B\ (x) . Since 
the variable y is cyclic under the Landau gauge, we can 
write the wave function as ip(x, y) = e lkyV ip(x), which in 
turn enables us to write (^) as 

Elpm = t 1pm~2 + Sm-l^m-1 + 8 m ^m+l + *oVWf2) 

(15) 

where 

5 m = i+e*"™ + t-e-* m (16) 

with 

Mm = Om + k y C. (17) 

Let us assume that a = p/q where p and q are coprimc 
integers and that T x (> 2) is an integer. Then, the period 

of 9m^ is q (2q) for even (odd) p and that of 9$ is 2T X 
if we set a = 1. Thus, the period of 9 m is given by 

f L.C.M.(g,2T x ) for even p 

\ L. CM. {2q,2T x ) for odd p ' 1 8j 

which leads to the relations [i m +T = Pm and S m +T = $m- 
Therefore, the Bloch condition along the x direction can 
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be expressed as ip m +T = e x ip m , and the characteristic 
matrix that gives rise to energy eigenvalues in presence 
of magnetic modulation can be written as 

\ 
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(19) 

where r\ = k x Tb. Noted that, when 7 = 0, the period of 
9 m , and thus the order of A, is given by To which equals 
to q (2q) for even (odd) p. 

We are now in position to discuss the formalism of 
Refs. 1 and 4. To this end, we refer to (14) and (16) 
of Ref. 1 as A F and ^ F , and (6) of Ref. 4 as A G , 
respectively. The most fundamental mistake made by 
the authors of Refs. 1 and 4 lies in that they em- 
ployed (@) instead of (Q) for developing the formalism. 
Since (0) is meaningful only in the % — » limit, it 
is evident that the formalism in Refs. 1 and 4 does 
not suit in describing the quantum-mechanical behav- 
ior of Bloch electrons in the presence of the magnetic 
field. To put it concretely, /x^ is an erroneous expres- 
sion as a result of the mistake; 9m ' (= 27rma) and 
6m F \= ~ (a^/(3) cos(27rm/3)] should be replaced by 9$ 

and 6m) , respectively, in order to describe pertinently the 
effect of the magnetic field given by (|l]) . Note that fi m 
in turn has a crucial influence on the order of the char- 
acteristic matrix. Since the period of 6m' F ' G ^ is q and 

F becomes T' = qT x , 
Thus, the order of A F 



( l'F) 

that of 9m is T x , the period of /i 



resulting in ip m +T' = c 



ik^T'b 



''Pi, 



is T" for generic values of q and T x . However, as shown 
above, the order of the corrected matrix A is not T" but 
T. Besides, the order of A under 7 = is not q but 
To, unlike to the assertion of Ref. 4. Another indication 
that the formalism of Refs. 1 and 4 is incorrect can be 
found in that A F and A are non-Hermitian, contrary 
to the general idea that a matrix which gives rise to real 
eigenvalues should be Hermitian. All these arguments in 
this paragraph indicate that the relevant numerical data, 
i.e., figures 3 and 4 in Ref. 1 and all the figures in Ref. 4, 
obtained by diagonalizing A F and/or A are not in the 
least the solutions of the problem that the authors of 
Refs. 1 and 4 attempted to solve. 

Figures 1(a) and 1(b) show plots of the energy eigen- 
values versus a in the absence of magnetic modulation 
(i.e., 7 = 0). We can see that figure 1(a) is the same as 
figure 3 of Ref. 6 while quite different from figure 1 of 
Ref. 4, which shows how terribly the energy spectrum is 
changed when 9m' F instead of 6^} is employed in the 
calculation. Indeed, even though the authors of Ref. 4 
were aware that their result differs from that of Ref. 6, 



they might have misunderstood that the two worksii 
deal with exactly the same problem. Comparison of fig- 
ure 1(b) with figure 3 of Ref. 4 also shows how sensitively 
the energy spectrum is affected by 6$ in the presence of 
the hopping anisotropy. Figures 2(a) and 2(b) show plots 
of the energy eigenvalues versus a in the presence of the 
magnetic modulation, where the same parameters as in 
Ref. 1 are chosen for comparison. The plots show that in- 
troducing magnetic modulation leads to very complicated 
subband structure; gap closing and subband broadening 
observed in the case of a square latticeBu can also be seen. 
Comparison of figures 2(a) and 2(b) with figures 3 and 4 
of Ref. 1 also shows the importance of the correct choice 
of 9m^ ■ For instance, the straight lines of E — —2 and/or 
E = —4 along the a axis in some plots of Refs. 1 and 
4 do not appear if we use correct form of the magnetic 
phase factor [i.e., (p))] in the calculation. 

Before ending up this comment, we would like to men- 
tion two remarks. One is that, since we focused our at- 
tention in this comment on the formalism of Refs. 1 and 
4, detailed description of the effect of magnetic modu- 
lation on the energy spectrum of the hexagonal lattice 
is still lacking and thus further study on this problem 
is required for better understanding. The other is that 
(O) and thus @ can also be obtained without difficulty 
if we starts directly from the tight-binding Hamiltonian 
given by 



(20) 



where 9. 
checked that 



A • dl. Indeed, it can be easily 



'm,n;m±2,n — ^m,n;m+l,n±l — i^m^ 

i — l,n±l = ±6 m -l, (21) 



y m.n;m- 



and that combining (21) with j20| ) leads to (|l3|), which 
confirms the correctness of the formalism derived in this 
comment. 
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